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Although greenhouse gases absorb primarily long-wave radiation, 

they also absorb short-wave radiation. Recent studies have highlighted 
theimportance of methane short-wave absorption, which enhances 

its stratospherically adjusted radiative forcing by up to - 15%. The 
corresponding climate impacts, however, have been only indirectly 
evaluated and thus remain largely unquantified. Here we present a 
systematic, unambiguous analysis using one model and separate 
simulations with and without methane short-wave absorption. We find 
that methane short-wave absorption counteracts -30% of the surface 
warming associated with its long-wave radiative effects. An even larger 
impact occurs for precipitation as methane short-wave absorption offsets 
~60% of the precipitation increase relative to its long-wave radiative 
effects. The methane short-wave-induced cooling is due largely to cloud 
rapid adjustments, including increased low-level clouds, which enhance 
the reflection of incoming short-wave radiation, and decreased high-level 
clouds, which enhance outgoing long-wave radiation. The cloud responses, 
inturn, are related to the profile of atmospheric solar heating and 
corresponding changes in temperature and relative humidity. Despite our 
findings, methane remains a potent contributor to global warming, and 
efforts to reduce methane emissions are vital for keeping global warming 
well below 2 °C above preindustrial values. 


® Check for updates 


The atmospheric concentration of methane (CH,) has increased 
by about a factor of 2.4 since preindustrial times (from -0.75 to 1.8 
parts per million by volume (ppm)), resulting in an effective radia- 
tive forcing (ERF; Methods) of 0.496 + 0.099 W m? (from 1850 
to 2019)!, with similar estimates based on the stratospherically 
adjusted radiative forcing (SARF; Methods)* *. Due to methane’s 


potency as a greenhouse gas (its global warming potential is 27.9 
times that of CO, on a 100 yr time horizon’), its relatively short 
lifetime (-1 decade) and chemical reactions in the atmosphere (for 
example, tropospheric ozone production), considerable interest 
exists in targeting CH, emissions to mitigate climate change and to 
improve air quality** . 
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Recent studies” have highlighted the importance of CH, 
short-wave (SW) absorption at near-infrared wavelengths—which is 
lacking in many climate models'—resulting in up to an -15% increase 
in its SARF compared with the long-wave (LW) SARF”. A more recent 
study? found a smaller increase in SARF, at 7%, which was attributed, in 
part, to the inclusion of CH, absorption of solar mid-infrared radiation 
in the 7.6-m-band spectral region. The reduced forcing is because 
thisspectral region impacts mainly stratospheric absorption. CH, SW 
absorption has regional ‘hotspots’, including near bright surfaces (for 
example, deserts) and above clouds (for example, oceanic stratus cloud 
decks)”. Such bright regions enhance the upward reflection of sunlight, 
which inturn enhances top-of-the-atmosphere (TOA) and tropopause 
CH,SW instantaneous radiative forcing (IRF). Considerable uncertainty 
exists, however, asthis forcing depends on several quantities, includ- 
ing the cloud radiative effect, CH, vertical profile” and surface albedo 
specification”. In particular, large spatial gradients in the SW forcing 
are caused by near-infrared surface albedo”. 

These studies focus largely on how CH, SW absorption impacts its 
radiative forcing, with some also addressing the corresponding rapid 
adjustments (ADJs, surface-temperature-independent responses). For 
example, CO, and CH, (with SW absorption) fixed sea surface tempera- 
tures (SST) and slab ocean simulations were compared to show that 
ADJs associated with CH, SW radiative effects act to mute precipitation 
increases“ due to enhanced warming of the upper troposphere and 
lower stratosphere (UTLS). Methane SW ADJs were also investigated 
in Precipitation Driver and Response Model Intercomparison Project 
(PDRMIP)” simulations. Models that lack CH, SW absorption yield a 
positive overall ADJ (acting to increase net energy into the climate sys- 
tem), whereas models that include CH, SW absorption yield a negative 
overall ADJ (acting to increase net energy out ofthe climate system)". 
This difference is due to a more negative tropospheric temperature 
adjustment and negative as opposed to positive stratospheric and 
cloud adjustments in models that include CH, SW absorption. These 
negative adjustments, inturn, are consistent with stronger UTLS warm- 
ing, which promotes enhanced outgoing LW radiation to space and 
high-level cloud reductions, which further promote enhanced outgoing 
LW radiation. Although other model differences (for example, cloud 
parameterizations and CH, vertical profile) may impact this result, 
the implication is that CH, SW absorption may not lead to additional 
surface warming. 

Although the importance of CH, SW absorption has been recog- 
nized, a comprehensive (and systematic) analysis of how it impacts 
the climate system remains to be conducted. Inthis study, we perform 
experiments to rigorously assess CH, SW radiative impacts on the 
climate system, including ADJs, surface-temperature-mediated feed- 
backs and the overall climate response. 


Results 

Asuite ofidealized methane-only time-slice perturbation simulations 
(Table 1 and Methods) are conducted with the National Center for 
Atmospheric Research Community Earth System Model version 2.1.3 
(CESM2)?. CESM2 includes the newest model components, including 
the Community Atmosphere Model version 6 (CAM6). Unlike many cli- 
mate models”, CAM6 includes CH, SW absorption in the near-infrared 
bands except the mid-infrared band in its radiative transfer parameteri- 
zation. For each methane perturbation (2x, 5x and 10x preindustrial 
atmospheric CH, concentrations) considered and ocean boundary 
condition (fixed climatological sea surface temperatures, fSST, ver- 
sus coupled ocean), we conduct pairs of identical experiments, one 
that includes CH,LW + SW radiative effects and one that lacks CH, SW 
radiative effects (Table 1). This allows quantification of the response 
signals (relative to preindustrial CH,) to CH, LW + SW, LW and SW radia- 
tive effects, abbreviated as CHa su, CH4,,, and CH, respectively. 
Under radiative transfer experiments (Methods), CH, radiative effects 
include the IRF only (the initial perturbation to the radiation balance). 


Table 1 | Description of CESM2/CAM6 methane experiments 


Experiment Description 
10xcHEXP 10xCH, with CH, LW+SW radiative effects? 
10xCHEXPsw 10xCH, with CH, SW radiative effects turned off 
(LW effects)? 
pie Preindustrial CH, with CH, LW+SW radiative effects 
pIEXP. Preindustrial CH, with SW radiative effects turned off 
NOSW 
(LW effects) 
Signal Description 
10xCH awasw = TOXCHEXP- PIP? Response to CH, LW+SW 


radiative effects” 


Response to CH, LW 


10xCHaiw = 10xCH iNosw7PINOsw Rc b 
radiative effects 


10*CHasw = 
(10xCHEXP- PI?) -(10xCH §Nosw7 PI Nosw) 


Response to CH, SW 
radiative effects” 


Top half of the table displays model experiments; bottom half displays response signals 
(difference of experiments). Atmosperic methane concentrations are perturbed relative to 
the preindustrial (Pl=year 1850) concentration; 10xCH, (preindustrial) represents 7.9 (0.79) 
ppm. Experiments are performed with both fixed climatological sea surface temperatures 
and a coupled ocean. The former allows quantification of the ADJs/fast responses; the latter 
allows quantification of the total climate response. The difference (total climate response 
minus fast response) quantifies the slow, surface-temperature-mediated feedback response. 
*Analogous sets of simulations are also conducted for both 5xCH, (3.95 ppm) and 2xCH, 
(1.58 ppm). ^Radiative effects include IRF for radiative transfer simulations; IRF and ADJs— 
which can be LW and SW adjustments—for fSST simulations; and IRF, ADJs and feedbacks for 
coupled simulations. 


Under fSST experiments, CH, radiative effects can induce an ERF, which 
includes both the IRF and ADJs (change in state in response to IRF, but 
excluding changes in sea surface temperatures). Rapid adjustments 
can be LW adjustments (for example, tropospheric and stratospheric 
temperatures), SW adjustments (surface albedo) or both SW and LW 
adjustments (clouds). The coupled ocean-atmosphere experiments 
quantify the total climate response, including the IRF, ADJs and the 
slow, surface-temperature-mediated effects. 


Methane SW versus LW total climate responses 
Figure 1 shows the global mean change in near-surface air tempera- 
ture and precipitation in coupled ocean-atmosphere CESM2 simu- 
lations, which is the total response (including IRF, adjustments and 
surface-temperature-mediated feedbacks) to increases in atmospheric 
methane concentrations, including 2xCH,, 5xCH, and 10xCH, relative 
to preindustrial (see also Extended Data Figs. 1 and 2). For all three 
perturbations, CH4,,—which represents the total climate response 
to methane LW IRF, adjustments and feedbacks-yields an increase 
in near-surface air temperature and precipitation (warming and 
‘wetting’). Significant global mean warming of 0.09, 0.68 and 1.24 K 
occurs for 2xCHgw, 5*CH4,w and 10xCH, y; similarly, global wetting 
of 0.001 (not significant at the 90% confidence level), 0.035 and 
0.063 mm d? occurs (corresponding precipitation changes are 0.04, 
1.2 and 2.1%). Interestingly, CH4,w.sw—which represents the total climate 
response to methane LW + SWIRF, adjustments and feedbacks—yields 
muted warming and wetting (except for 2xCH4,w.sw). This is due to 
SW effects (including the IRF, adjustments and feedbacks), where 
significant global cooling occurs for 5xCHysy and 10xCH,s, at -0.23 
and -0.39 K. Similarly, a significant decrease in global mean precipi- 
tation occurs under these two methane perturbations at -0.021 and 
-0.039 mm d” (-0.7 and -1.3%). Most of the precipitation decrease 
occurs over tropical oceans (for example, Extended Data Fig. 2c). 
The decrease in precipitation is consistent with atmospheric ener- 
getic constraints—in the global mean, the primary balance is between 
net atmospheric radiative cooling and condensational heating from 
precipitation", As atmospheric SW absorption increases, net 
radiative cooling decreases, which is consistent with a decrease in 
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Fig. 1| Total climate response to methane perturbations. a,b, Global annual 
mean near-surface air temperature (a) and precipitation (b) response for 
2xCH, (0.79 ppm), 5xCH, (3.16 ppm) and 10xCH, (7.11 ppm) from coupled 
simulations (which include the IRF, adjustments and feedbacks). Responses are 
decomposed into CH, ws, CH4,,, and CH;sy. Also included are the least-squares 
regression lines (dotted). Solid circles represent a significant response at the 
90% confidence level, based ona standard ttest. The thin black vertical line 
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shows the present-day CH, perturbation of 1.1 ppm. c-e, The estimated (from the 
regressions) present-day CH, climate responses for near-surface air temperature 
(c), precipitation (d) and apparent hydrological sensitivity (e). Except for 

e, the CH, u.s, bar is equal to the sum ofthe CH,,,, and CH,sw bars. Error bars 
inc-eshow the1 s.d. uncertainty estimate of the regression slope, which is 
estimated from the three like-coloured data points (CAM6 methane simulations) 
in aand b. The error-bar centre is the regression-estimated response. 


precipitation. Except for the 2xCH, perturbation, CH;s, offsets over 
~30% of the surface warming and -60% of the wetting associated with 
CH, ; that is, SW absorption offsets twice as much ofthe precipitation 
increase, compared with the surface warming. 

We estimate the present-day CH, climate response (ACH, of 
1.1 ppm) from least-squares regressions applied to our idealized 
2xCH,, 5xCH, and 10xCH, simulations (Fig. 1a,b). Figure 1c shows 
the corresponding near-surface air temperature response, decom- 
posed into CH, w.sy CHapw and CHysy. We find global warming of 
0.17 K in response to present-day CH, (relative to preindustrial); this 
is decomposed into warming of 0.20 K from CH,,,, and -0.04 K (cool- 
ing) from CH,sy. Our estimate of 0.17 K for CH4,,.swis less than that 
given in the newest IPCC report (based on 2019 relative to 1750 and a 
two-layer emulator) at 0.28 K, with a 5-95% range of 0.19 to 0.39 K (ref. ') 
(discussed in Supplementary Note 1). For global mean precipitation 
(Fig. 1d), precipitation increases of 0.16% and 0.31% occur for CH, w.sw 
and CH, respectively; a precipitation decrease of -0.15% occurs 
under CH,sw. The apparent hydrological sensitivities (defined as the 
changein precipitation divided by the change in surface temperature)” 
are 0.97, 1.51 and 4.30% K“ for CH4, usw, CH4rw and CHysw, respectively 


(Fig. le; discussed in Supplementary Note 2). Our decomposition helps 
to explain the larger apparent hydrological sensitivity to methane 
foundin PDRMIP models (many of which lack CH, SW radiative effects; 
discussed in Supplementary Note 2). 


Radiative flux components 

To understand the cause of the CH,4, surface cooling in coupled simula- 
tions, we evaluate the radiative flux components (Methods)- including 
ERF, IRF and the ADJs—in the fSST experiments. Figure 2a shows the 
TOA radiative flux components in response to 10XCH, w.sw, LOxCHaw 
and10xCH,s,. The IRF is 2.08 W m”, with 1OxCH4,y and 10XCH,s, both 
contributing positive values at 1.81 and 0.27 W m”, respectively. Thus, 
the 10xCH,sw IRF increases the 10xCH,,, IRF by 15% (13% for 5xCH, 
and 2xCH,). A previous study found a similar 15% increase under a 
750-1,800 ppb CH, perturbation". A smaller increase of 696 was found 
at the tropopause”, but the partitioning of SW IRF and LW IRF at the 
tropopause will differ from the TOA”. We also note that the presence of 
clouds increases the 10xCH,sw IRF from 0.20 W m under clear-sky con- 
ditions to 0.27 W m? under all-sky conditions (a 35% increase; 5xCH,s, 
and 2xCHysy yield 27% and 33% increases, respectively). The increased 
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Fig. 2| TOA radiative flux components and ADJ decomposition for 10xCH,. 
a, Global annual mean TOA ERF, IRF and ADJ. b, Global annual mean TOA surface 
temperature, tropospheric temperature, stratospheric temperature, water 
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vapour, surface albedo, cloud and total ADJ for 10xCH,. Responses that are 
not significant, based on a standard t test at the 90% confidence level, have 
unfilled bars. 


forcing due to clouds is related to increased absorption path lengths 
in the CH, bands caused by multiple scattering®”. 

The 10xCH,sw4 y and 10xCH,,y ERFs are also positive, at 1.69 
and 2.13 W m”, respectively, but the 10xCH,sw ERF is negative, at 
-0.44 W m”. Thus, 10xCH,,, acts to reduce the 10xCH,,, ERF by 
21%. The difference between ERF and IRF is due to ADJs; 10XCH,swaw 
yields a negative ADJ, at -0.40 W m”, which is due to the 10xCH,.y of 
-0.77 W m? (relative to the positive ADJ for I10xCH, of 0.37 W m?). 
Thus, 10xCH,s, drives a strong negative ADJ, offsetting its smaller posi- 
tive IRF (by a factor of -3), leading to a negative ERF of -0.44 W m”. This 
10xCH,sy negative ERF is consistent with the corresponding decrease 
in near-surface air temperature previously discussed (Fig. 1). We note 
that some ofthe 10xCH,sw adjustments are LW adjustments (discussed 
in ‘ADJ decomposition’). 

Qualitatively similar results are obtained from 5xCH, and 2xCH, 
(Supplementary Note 3 and Extended Data Fig. 3). Atmospheric and 
surface radiation contributions to the TOA radiation (ERF) changes 
are discussed in Supplementary Note 4 (see also Supplementary Fig. 1). 


ADJ decomposition 

To further understand the ADJs and climate impacts of CH, SW absorp- 
tion, Fig. 2b shows the decomposition of TOA ADJs (for 10xCH,) into 
the tropospheric temperature, stratospheric temperature, surface 
temperature, water vapour, albedo and cloud adjustment (Methods). 
Clouds are the main driver of the relatively large negative 10xCH,sw ADJ. 
The corresponding cloud adjustment is -0.58 W m, which is 75% of the 
total ADJ. The stratospheric temperature adjustment also contributes, 
at -0.15 W m”, as does the tropospheric temperature adjustment, 
at -0.11 W m”. The water-vapour adjustment—at 0.10 W m ?-acts to 
oppose these negative adjustments. The remaining ADJs, including 
surface temperature and albedo, are relatively small. Similar results 
are obtained for 5xCH, and 2xCH, (Extended Data Fig. 4). 

The 10xCH,s, cloud adjustment (Fig. 2b) is due to both SW radi- 
ation, at -0.42 W m” (Extended Data Fig. 4c), and LW radiation, at 
-0.16 W m” (Extended Data Fig. 4b). The corresponding 10xCHysw 
temperature and water-vapour adjustments are consistent with 
atmospheric warming (particularly in the UTLS; Fig. 3b), which leads to 
enhanced outgoing LW radiation (a negative LW adjustment; Extended 
Data Fig. 4b); the warming likewise increases water vapour (a green- 
house gas), which acts to decrease outgoing LW radiation (a positive LW 


adjustment; Extended Data Fig. 4b). Supplementary Note 5 discusses 
the decomposition of surface (and atmospheric) ADJs for 10xCH, 
(see also Supplementary Fig. 2). 

Recent analyses? have shown similar results across different 
kernels”, including the Geophysical Fluid Dynamics Laboratory 
kernel used here. Nonetheless, we repeat our ADJ calculations with 
the CloudSat/CALIPSO” radiative kernel and find similar results 
(Supplementary Note 6 and Extended Data Fig. 5). 


Understanding the cloud adjustment 

The negative 10xCH,sw cloud adjustment- including negative TOA 
SW and LW contributions—is consistent with the change in the global 
mean vertical profile of cloud cover (Fig. 3d, dashed line). This includes 
increased low-level cloud cover (peaking at 800 hPa) and enhanced 
reflection of SW radiation (a negative adjustment) but decreased 
high-level cloud cover (peaking at 100 hPa) and enhanced outgoing 
LW radiation (a negative adjustment). The change in the vertical profile 
of cloud cover is related to the change in relative humidity (RH), which 
increases below -500 hPa but decreases aloft (Fig. 3c, dashed line). The 
corresponding correlation, r, from the surface up to the lower strato- 
sphere (up to -100 hPa) is 0.86, suggesting an increase (decrease) in 
RHisassociated with more (fewer) clouds. The change in RH is consist- 
ent with the change in the vertical profile of temperature (r = -0.76; 
Fig. 3b, dashed line), which in turn is related to the atmospheric SW 
heating rate (r = 0.88; 3a dashed line). Thus, we suggest the 10xCHysw 
cloud responseis driven ultimately by the atmospheric SW heating-rate 
profile, which decreases in the low/mid troposphere (below -700 hPa) 
but increases aloft, peaking in the UTLS at 100 hPa. This is consistent 
with the traditionally defined aerosol-cloud semi-direct effect” ^*, 
whereby solar heating (for example, from black carbon) increases 
atmospherictemperature and decreases RH, leadingto cloud burn-off 
(with the opposite occurring in the lower troposphere). Atmospheric 
cooling below -800 hPa and warming aloft also imply an increase in 
stability, which is also probably associated with the increase in low 
cloud cover. Similar responses occur under 5xCH, (Supplementary 
Fig.3) and (although weaker) 2xCH, (Supplementary Fig. 4). 


Atmospheric SW heating response profile 
Theglobal annual mean CH, instantaneous SW heating-rate response 
profile is not related to the vertical profile of the CH, concentration, 
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Fig. 3 | Global annual mean vertical profiles of fast responses for 10xCH,. 
a-d, Atmospheric SW heating rate (QRS) (a), temperature (T) (b), relative 
humidity (RH) (c) and cloud fraction (CLOUD) (d) for 10xCH,. Panels include the 
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contributions from 10XCH,,w.sw, L1OXCHaw and 10xCH;sw. Solid dots represent 
asignificant response at the 90% confidence level, based on a standard t test. 
Also included ina is the instantaneous SW heating-rate profile (1OXCH,sy ir). 


which in CESM2/CAMG has a uniform distribution in the troposphere 
(up to ~200 hPa) and then exponentially decreases aloft (Extended Data 
Fig. 6), consistent with chemical destruction of CH, above the tropo- 
pause. Instead, the instantaneous SW heating-rate response profile is 
related to overlap of the three CH, SW absorption bands with water 
vapour. Under clear-sky conditions, with water-vapour SW absorption 
inthethree methane SW bands (Methods) turned off (using the Parallel 
Offline Radiative Transfer (PORT) model), the vertical profile of CH, 
SW instantaneous absorptionis relatively uniform in the troposphere, 
peaking in the UTLS (Fig. 4a). Adding back the SW absorption by water 
vapour leads to the characteristic SW heating-rate response profile 
(as in Fig. 3a), with decreases in the lower troposphere and increases 
aloft, peaking in the UTLS. As expected, the 10xCH,.y clear-sky IRF 
increases (from 0.20 to 0.40 W m?) when the overlapping SW absorp- 
tion by water vapour is turned off. 

Since water vapour is at its maximum in the lower troposphere, 
these SW absorption bands are already highly saturated in the lower 
atmosphere at preindustrial CH, concentrations, so perturbing meth- 
ane does not lead to an increase in SW heating here. However, methane 
SW radiative effects enhance SW absorption aloft (increase in SW 
heating rate). This reduces the amount of solar radiation in these three 
bands that can be subsequently absorbed by water vapour inthe lower 
troposphere, which results in the SW heating-rate decrease below 
~700 hPa. Similar results are obtained under all-sky conditions (Fig. 4b). 
The 10xCH,sw IRF increases (from 0.27 to 0.43 W m?) when the overlap- 
ping SW absorption by water vapour is turned off. Here, however, even 


with water-vapour SW absorption (in the three methane bands) tuned 
off, there is still a decrease in the instantaneous SW heating rate near 
800 hPa. This appears to be related to clouds, which peak at about the 
same level. Extended Data Fig. 7 shows similar plots but based on three 
different latitude bands: the low latitudes (30° S-30? N); mid-latitudes 
(30° N-60° Nand 30° S-60° S) and the high latitudes (60° N-90° Nand 
60° S-90? S). There are some differences relative to the global mean 
(Fig. 4), but the results are generally similar. For example, absorption by 
water vapour (in the three methane bands) is more important in the low 
latitudes (Extended Data Fig. 7e), consistent with the larger amount of 
water vapour (specific humidity) in the tropics. To summarize, methane 
SW instantaneous radiative effects result in a vertical redistribution of 
atmospheric SW heating, with enhanced SW heating aloft (maximizing 
inthe UTLS) but decreased SW heating in the lower troposphere. This, 
in turn, leads to the corresponding cloud-cover changes (increased 
low-level but decreased high-level cloud cover) and negative cloud 
adjustment, through modification of atmospheric temperature and 
relative humidity. 


Climate feedbacks under methane SW radiative effects 

Figure 5 shows the radiative kernel decomposition applied to the 
coupled ocean-atmosphere simulations for 10xCH,s, (5*CHysy and 
2xCH,sy are included in Supplementary Fig. 5). We also include the 
previously discussed ADJs (‘fast’ responses from the fSST runs) and 
the difference between the coupled and fSST decompositions (the 
surface-temperature-induced 'slow' feedbacks). Note that we do not 
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near-infrared bands that methane absorbs in. Also includes the climatological 
cloud fraction (CLOUD). 
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Fig. 5| TOA radiative flux decomposition for the total response, ADJ and 
feedback for 10xCH,sw. Global annual mean TOA surface temperature, 
tropospheric temperature, stratospheric temperature, water vapour, surface 
albedo, cloud and total radiative flux decomposition for 10xCH,sw. The first bar 
in each like-coloured set of three bars represents the total response (from the 
coupled ocean simulations), the second bar represents the ADJ (fast response) 
andthethird bar represents the surface-temperature-induced feedback (slow 
response). Responses not significant, based ona standard ttest at the 9096 
confidence level, have unfilled bars. Units are W m. 


normalize our feedbacks by the change in global mean surface tem- 
perature; unnormalized feedbacks facilitate comparison with the 
ADJs. Thus, positive/negative feedbacks have the same meaning as 
positive/negative ADJs (positive is an increase in net energy; negative 
is a decrease in net energy). 

In most cases, the slow feedback dominates the sign of the overall 
response, consistent with the climate system acting to restore TOA 
radiative equilibrium. For example, the slow tropospheric temperature 
feedback is positive at 1.14 W m (which is offset to some extent by 
the water-vapour feedback at -0.62 W m ?). Both of these feedbacks 
are consistent with tropospheric cooling (Supplementary Fig. 6b). 


For clouds, however, the ADJ and the slow feedback are both negative, 
with a larger value for the ADJ, at -0.58 versus -0.37 W m”. Thus, the 
surface cooling in response to CH, SW radiative effects is due largely 
to cloud ADJs, but surface-temperature-induced cloud feedbacks also 
act to cool the planet. 

The 10xCH,s, cloud feedback is dominated by increases in 
low-level (and mid-level) clouds, with weaker decreases in high-level 
clouds (Supplementary Fig. 6d and Supplementary Note 7). Similar 
results exist for 5xCH, wy (Supplementary Fig. 7), but weaker results 
exist for 2xCH,sw (Supplementary Fig. 8 and Supplementary Note 7). 


Conclusions 

Using targeted climate model simulations, we have shown that methane 
SW absorption and the associated ADJs act to reduce its ERF by ~20% 
and mute its warming and wetting effects in coupled simulations by 
up to 30% and 60%, respectively. Similar simulations with additional 
climate models are needed to understand the robustness ofthe results 
presented here- particularly since the CH, SW IRF is dependent on 
uncertain quantities, such as the cloud radiative effect”, the surface 
albedo??? and the CH, vertical profile”. However, the indirect assess- 
ment of multiple models from PDRMIP CH, simulations supports our 
findings”. In fact, expanding on the results of ref. ?^, we find a 2096 
decrease in ERF, 45% less warming and 65% less wetting in models that 
include CH, SW absorption versus those that do not (Supplementary 
Note 8 and Extended Data Fig. 8). 

Althoughthe SW radiative effects associated with the present-day 
methane perturbation remain relatively small, they could be quite 
large by the end of the century—Shared Socioeconomic Pathway 3-7.0, 
which lacks climate policy and has ‘weak’ levels of air-quality control 
measures”, features end-of-century increases of CH, concentrations 
approaching 5x preindustrial (3.4 ppm). Overall, methane remains a 
potent contributor to global warming, and emissions reductions are 
a vital component of climate change mitigation policies and for con- 
tinued pursuit ofthe climate goals laid out under the Paris Agreement. 


Online content 

Any methods, additional references, Nature Portfolio reporting sum- 
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Methods 

Radiative-forcing definitions 

The IRF is the initial perturbation to Earth's radiation budget and does 
not account for ADJs. We diagnose IRF using the PORT model", which 
isolates the Rapid Radiative Transfer Model for general circulation 
models (RRTMG)**” radiative transfer computation from the CESM2- 
CAM6 model configuration (more details on RRTMG are presented in 
‘CESM2/CAM6 simulations’). PORT simulations are run for 16 months; 
the last 12 months are used to diagnose annual mean IRF. PORT is also 
usedto verify our methodology to remove RRTMG CH, SW absorption 
(the SW IRF is zero in the CHEX?.., and PIE, experiments, and the LW 
IRF isunchanged). 

TheERF is definedasthe net TOA radiative flux difference between 
the perturbed and base simulation, with climatological fixed SSTs and 
sea-ice distributions and no correction for land surface-temperature 
change?" We note that the contribution of land surface warming/ 
cooling to the ERF in our simulations is relatively small (<5% of the 
ERF; Supplementary Note 9). ERF can be decomposed into the sum of 
IRF and the ADJs. 

The SARF is equal to the sum of the IRF and the stratospheric 
temperature adjustment. Thus, the difference between ERF and SARF 
is that ERF includes all adjustments, whereas SARF includes only the 
adjustment due to stratospheric temperature change^*"?, 


CESM2/CAM6 simulations 

We conduct pairs of identical simulations, one that includes CH, 
LW+SW radiative effects (CH§**) and one that lacks CH, SW radiative 
effects (CHEX? sw; Table 1). The latter simulations are conducted by 
turning off CH, SW absorption in the three near-infrared bands-1.6- 
1.9 um, 2.15-2.50 um and 3.10-3.85 um- in CAM6’s radiative transfer 
parameterization (RRTMG). RRTMG does not include methane SW 
absorption in the mid-infrared band at 7.6 pm. The sign ofthe CH, SW 
IRF (at the tropopause) depends on the increased absorption in the 
troposphere since the downward SW flux at the tropopause is always 
decreased due to absorption in the stratosphere”. Including the 7.6 jum 
band primarily increases CH, SW absorption in the stratosphere”. This 
reduces the forcing from the downward irradiance, with negligible 
change to the forcing from the upward irradiance; that is, the tropo- 
pause SW IRF is reduced. Thus, if RRTMG included the 7.6 um methane 
band, we would expect the CH, SW IRF at the TOA to increase due to 
the increase in stratospheric absorption. This, however, will result in 
alarger (negative) stratospheric temperature adjustment. 

RRTMGisan accelerated and modified version of RRTM and uses 
the correlated k-distribution method to treat gas absorption”. RRTMG 
calculates irradiance and heating rate in broad spectral intervals 
while retaining a high level of accuracy relative to measurements and 
high-resolution line-by-line models. Sub-grid cloud characterization 
is treated in both the LW and SW spectral regions with the Monte Carlo 
Independent Column Approximation” using the maximum-random 
cloud overlap assumption. RRTMG divides the solar spectrum into 14 
SW bands that extend over the spectral range from 0.2 pm to 12.2 pm. 
The infrared spectrum in RRTMG is divided into 16 LW bands that 
extend over the spectral range from 3.1 jum to 1000.0 um. 

Few studies have evaluated broadband radiative transfer codes 
against benchmark calculations, particularly for CH, SW IRF. This is 
in part because the radiation parameterization in many climate mod- 
els lacks an explicit treatment of CH, SW absorption™"®. The 6-band 
SOCRATES SW spectral file configuration used in the Met Office Unified 
Model significantly underestimates CH, SW tropopause and surface IRF 
by around 45% compared with the 260-band configuration”. Similarly, 
RRTMG-the radiative transfer model used here—was recently found 
to underestimate CH, (and CO;) SW IRF by 25-45% (ref. **). This implies 
that there are opportunities for improvement in the parts of the spec- 
trum where the absorption by these gases is weak but not zero. Thus, 
incorporating CH, SW absorption in more models’ radiative transfer 


codes is only part ofthe solution—making sure their radiative transfer 
codes have a validated treatment of SW absorption by CH, (and other 
greenhouse gases) is also vital. We also note that N,O is not represented 
inthe SW part of RRTMG. 

The Community Land Model version 5? provides both the surface 
albedo, area-averaged for each atmospheric column, and the upward 
LW surface flux, which incorporates the surface emissivity, for input 
to the radiation. For the SW, the surface albedos are specified at every 
grid point at every time step. The albedos are partitioned into two 
wavebands (0.2-0.7 um and 0.7-12.0 pum) for both direct and diffuse 
incident radiation 6. Surface albedos for ocean surfaces, geographi- 
cally varying land surfaces and sea-ice surfaces are distinguished. They 
depend onthe solar zenith angle, the amount and optical properties of 
vegetation and the optical properties of snow and soil”. 

Rapid adjustments—which can be SW or LW adjustments—are 
estimated by subtracting the preindustrial control (PI^) fSST experi- 
ment from each perturbation fSST experiment. For example, to quan- 
tify the ADJs in response to a tenfold increase in preindustrial 
atmospheric methane concentration, we take the 10xCH, fSST simula- 
tion minus the preindustrial fSST simulation (IOxCH7** - PIP®). This 
signal (IOxCH,, ws) includes the methane LW+SW IRF and its impact 
on LW and SW adjustments under the fSST boundary condition 
(ERF = IRF + adjustments). ADJs due to CH, LW IRF and its impact on 
LW and SW adjustments (10xCH,,,) are estimated from 10xCH EXP sw- 
PIEX w. Similarly, ADJs due to CH, SW IRF and the impact of CH, SW 
absorption on LW and SW adjustments (10xCH,sw) are estimated from 
(10xCHẸ® - PI^) — (10xCH RR sw- PIER). Specific details on how the 
ADJs are estimated (via radiative kernels) are discussed in ‘Calculation 
of ADJs’. A similar procedure is used to quantify the total climate 
impacts from the coupled ocean simulations, which include the IRF, 
adjustments and surface-temperature-mediated feedbacks. 

We note that an alternative experimental design where methane 
LW radiative effects are removed could be implemented. As our goal is 
to understand the impacts of adding CH, SW absorption (which many 
models lack) to the LW forcing (which models already have), our experi- 
mental design is based on the all-but-one type of experimental design. 
Our simulations therefore target the inclusion of CH, SW absorption, 
allowing quantification of its associated ADJs and climate impacts. This 
isin contrast to the studies discussed, which evaluate CH, SW radiative 
effects either by contrasting CH, versus CO, (which lacks strong SW 
absorption) simulations!'é or by comparing models that include CH, 
SW absorption with models that do not”. In the latter, other model 
differences (for example, cloud parameterizations and CH, vertical 
profile) may beimportant. 

All CESM2/CAM6 simulations are conducted with a 1.9? x 2.5? 
horizontal resolution and 32 vertical levels in the atmosphere. Fixed 
SST experiments are run for 32 years each, the last 30 of which are 
used to quantify ERF and the ADJs/fast responses. Coupled ocean 
simulations are run for 90 years each, starting from a pre-spun-up 
preindustrial control simulation in year 321. The last 40 years of the 
coupled experiments—when the net TOA radiative flux stabilizes—are 
used to quantify climate impacts. The surface-temperature-mediated 
slow response is calculated as the difference between coupled ocean 
and fSST experiments". A 90-year coupled ocean simulation has not 
yet reached equilibrium, so we refer to these simulations as being in 
near equilibrium (computational cost restrictions prohibited longer 
integrations), similar to previous projects including PDRMIP?^5, Our 
CESM2/CAM6 simulations do not include interactive chemistry; we 
therefore do not address possible atmospheric chemistry implications 
(for example, changes in ozone and stratospheric water vapour) or 
changes in methane lifetime. 


Calculation of ADJs 
The ADJs (for example, clouds, water vapour and temperature) in 
the climatological fixed SST experiments are estimated using the 
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radiative kernel method*"*”*°, Radiative kernels represent the radia- 
tiveimpacts from small perturbations ina state variable (for example, 
temperature, water vapour and surface albedo). Subsequently, ADJs 
can be computed by multiplication of the kernel with the response of 
the state variable. We use the Python-based radiative kernel toolkit 
(downloaded from ref. ^) and the Geophysical Fluid Dynamics Labora- 
tory radiative kernel. 

Weusethe same radiative kernel procedureto calculate the unnor- 
malized (we do not divideby the change in global mean surface temper- 
ature) feedbacks. Specifically, surface-temperature-induced feedbacks 
are estimated by subtracting the ADJs (from fixed SST experiments) 
fromthe corresponding radiative kernel decomposition applied to the 
coupled experiments. Unnormalized feedbacks facilitate comparison 
with the ADJs. 

ERF can be decomposed as ERF = IRF + ADJ44 + ADJs; + ADJz5 + 
ADJwy + ADJ, + ADJ. + e, where ADJ,,; is the tropospheric temperature 
adjustment, ADJ,, is the stratospheric temperature adjustment, ADJ-s 
is the surface-temperature adjustment, ADJ,, is the water-vapour 
adjustment, ADJ, is the albedo adjustment, ADJ. is the cloud adjust- 
ment and e isthe radiative kernel error. Individual ADJs are estimated 

as ADJ, = T dx, where = is the radiative kernel and dxis the response 
of state variable x as simulated by CESM2/CAM6. Kernels are 
four-dimensional (latitude, longitude, pressure and month) fields 
for atmospheric temperature and specific humidity and three- 
dimensional (latitude, longitude and month) fields for surface tem- 
perature and surface albedo. Two sets of kernels are used: clear-sky 
kernels, where the fluxes are calculated without clouds, and all- 
sky kernels. 

As the radiative effect of clouds depends on several variables 
(fraction, ice and liquid-water content, droplet effective radius and 
so on), several approaches have been used to estimate cloud adjust- 
ments!*°°. Here we estimate cloud adjustments using the kernel 
difference method?, which involves a cloud-masking correction 
of cloud radiative-forcing diagnostics using the kernel-derived 
non-cloud adjustments and IRF according to ADJ, = (ERF - ERF,,) — 
(IRF - IRF.) — È x-tr. Ts, wv a ( ADJ; - ADJ xcs), where subscript ‘cs’ refers to 
clear-sky quantities. Thus, the kernel difference method relies onthe 
difference of all-sky and clear-sky kernel decompositions. See ref. ^ 
for additional details. 

The total ADJ is estimated as the sum of individual ADJs from the 
radiative kernel decomposition. Since we estimate IRF using PORT 
for all of our methane simulations, this can be used to estimate the 
radiative kernel error (e) ase = ERF - IRF — )..-17-1s, w,a,q(ADJ,). For exam- 
ple, the I0xCH,,,,, 4, ERF and IRF are 1.69 and 2.08 W m”, yielding an 
ERF - IRF difference of -0.39 W m”. Thesum of the individual ADJs from 
the kernel decomposition is -0.40 W m”. Thus, the radiative kernel 
error for IOXCH,,,, suis 0.01 W m”. Similar results hold for SXCH,,w.sw 
and 2xCH;ıw+sw, where eis 0.03 W m? and -0.02 W m”, respectively. 
Relative to the corresponding ERFs, these errors are «1.096, 3.1% and 
5.7%, respectively. As ref. ^ lacked an estimate of the IRF (which we 
estimate using PORT), they estimated € under select situations (where 
the SW or LW IRF is known to be zero). In these situations, they found 
that the residual term is small, being “6%, 12% and 2% of the ERF for 
10xBC LW, 3xCH, SW and 2%Solar LW in magnitude, respectively. The 
larger multimodel residual in the 3xCH, SW case is biased by a large 
relative residual in the HadGEM2 model, whereas residuals in the other 
four models analysed are close to 0.” Thus, our radiative kernel errors 
are relatively small, and comparable to those estimated from select 
PDRMIP simulations”. 

We note that methane IRF has an approximate square root 
dependency on concentration**!. PDRMIP 3xCH, simulations yield 
a 3xCH, IRF of 1.1 + 0.24 W m”, but nearly all of the PDRMIP models 
used year 2000 as the base year. This perturbation is thus similar to 
5-6x preindustrial CH, (our 5xCH, IRF is 1.18 W m”, with 0.14 W m? 
due to SW radiative effects). 


Statistical significance 

Statistical significance of a climate response is calculated using a 
two-tailed pooled ttest. An annual mean time series is calculated for 
both the perturbation experiment and the preindustrial base experi- 
ment (for example, at individual grid boxes or averaged globally), and 
their difference is taken. The null hypothesis of a zero difference is 
evaluated, with n; + n; - 2 degrees of freedom, where n; and n; are the 
number of years in the perturbation experiment and base (30 years 
forfSST experiments; 40 years for coupled ocean experiments). Here, 


(1; -1)52 (n -1)52 


the pooled variance, 2, is used, where S, and S, are the 


njtn-2 
sample variances. 

Asimilar procedure is used to quantify statistical significance of 
the radiative flux perturbations and rapid adjustments (for example, 
Fig. 2). These uncertainties are therefore relative to interannual vari- 
ability and do notaccountfor possible intermodel or kernel uncertain- 
ties (asin ref. ?, using 10+ PDRMIP models). As we have only one year of 
data for the IRF, we evaluate its uncertainty relative to the preindustrial 
base experiment with fixed SSTs. Nearly all of our ADJs under 10xCH, 
are significant at the 90% confidence level (the lone exception is the 
surface-temperature adjustment under 10xCH,.w). Similar conclusions 
also hold for 5xCH, (Extended Data Fig. 4d). Under 2xCH,, however, 
most of the ADJs under 2xCH,,, are not significant (Extended Data 
Fig. 4g), including the total ADJ. This is consistent with the relatively 
small 2xCH, SW IRF of 0.04 W m”. 

We also find similar results using an alternative kernel (CloudSat/ 
CALIPSO; Extended Data Fig. 5), so our ADJ conclusions are robust across 
these two kernels. Finally, we note that the ADJs in PDRMIP models that 
include CH, SW absorption (under 3xCH,, whichis a perturbation simi- 
lar to our 5x preindustrial CH,) are all significant at the 95% confidence 
level’®, and this includes the intermodel and kernel uncertainty. 


Data availability 

PDRMIP simulations can be accessed at https://cicero.oslo.no/en/ 
PDRMIP/PDRMIP-data-access. A core set of model data from our ideal- 
ized methane CESM2 simulations can be downloaded from Zenodo at 
https://doi.org/10.5281/zenodo.7596623. 


Code availability 
The Python-based radiative kernel toolkit and the GFDL radiative 
kernel can be downloaded from https://climate.rsmas.miami.edu/ 
data/radiative-kernels/. 
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Extended Data Fig. 1| Global maps of the total near-surface air temperature represent a significant response at the 90% confidence level, based on a standard 
response to methane perturbations. Annual mean near-surface air t-test. Climate responses are estimated from the coupled simulations, which 
temperature response for (a-c) IOxCH,, (d-f) 5xCH, and (g-i) 2XCH, decomposed include the IRF, adjustments and feedbacks. 


into (a, d, g) CHaiwssw, (b, e, h) CHypw and (c, f, i) CH,s4. Units are K. Solid dots 
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Extended Data Fig. 2 | Global maps of the total precipitation response to at the 90% confidence level, based ona standard t-test. Climate responses are 
methane perturbations. Annual mean precipitation response for (a-c) 10xCH,, estimated from the coupled simulations, which include the IRF, adjustments 
(d-f) 5xCH, and (g-i) 2xCH, decomposed into (a, d, g) CHaiwisw, (b, e, h) CHa and feedbacks. 


and (c, f, i) CHysy. Units are mm day”. Solid dots represent a significant response 
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Extended Data Fig. 3 | Top-of-the-atmosphere radiative flux components 

in response to methane perturbations. Global annual mean top-of-the- 
atmosphere (TOA) effective radiative forcing (ERF; black), instantaneous 
radiative forcing (IRF; green) and rapid adjustment (ADJ; blue) for (a-c) 10XCH,, 
(d-f) 5xCH, and (g-i) 2xCH, for (a, d, g) LW+SW, (b, e, h) LW and (c, f, i) SW 
radiative fluxes. The first bar in each like-colored set of three bars represents the 
contribution from CH, w;sw; the second bar represents CH4,w; and the third bar 


represents CH,s,. Responses not significant, based ona standard t-test at the 90% 
confidence level, have unfilled bars. These uncertainties are therefore relative to 
interannual variability. Units are W m”. ERF is estimated as the net TOA energy 
imbalance from the fixed SST experiments. IRF is estimated using an offline 
radiative transfer model, PORT. The total rapid adjustment is estimated as the 
sum of individual rapid adjustments from the radiative kernel decomposition. 
Panel (a) is identical to Fig. 2a from the main paper. 
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Extended DataFig. 4 | Top-of-the-atmosphere rapid adjustment 
decomposition in response to methane perturbations. Global annual mean 
top-of-the-atmosphere (TOA) surface temperature (purple), tropospheric 
temperature (cyan), stratospheric tem- perature (yellow), water vapor (red), 
surface albedo (orange), cloud (pink) and total (blue) rapid adjustment for (a-c) 
10xCH,, (d-f) 5xCH, and (g-i) 2xCH, for (a, d, g) LW+SW, (b, e, h) LW and (c, f, i) SW 
radiative fluxes. The first bar in each like-colored set of three bars represents the 


contribution from CH, sw; the second bar represents CH,,; and the third bar 
represents CHysw. Responses not significant, based on a standard t-test at the 90% 
confidence level, have unfilled bars. These uncertainties are therefore relative 

to inter- annual variability. Rapid adjustments are estimated from the radiative 
kernel decomposition applied to the fixed SST experiments. Units are W m”. 
Panel (a) is identical to Fig. 2b from the main paper. 
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Extended Data Fig. 5 | Top-of-the-atmosphere rapid adjustment 
decomposition using CloudSat/CALIPSO radiative kernels. Global annual 
mean top-of-the- atmosphere (TOA) surface temperature (purple), tropospheric 
temperature (cyan), stratospheric temperature (yellow), water vapor (red), 
surface albedo (orange), cloud (pink) and total (blue) rapid adjustment for (a) 
10xCH,, (b) 5xCH,, and (c) 2xCH,. The first bar in each like-colored set of three 


bars represents the contribution from CH, w.sw; the second bar represents 

CHA w; and the third bar represents CH;s,. Responses not significant, based 
onastandard t-test at the 90% confidence level, have unfilled bars. These 
uncertainties are therefore relative to interannual variability. Units are W m”. The 
total rapid adjustment is estimated as the sum of individual rapid adjustments 
fromthe radiative kernel decomposition. 
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Extended Data Fig. 6 | Global annual mean vertical profiles of the atmospheric methane concentration perturbation. Atmospheric CH, concentration [ppm] for 


10xCH, (black), 5xCH, (blue) and 2xCH, (green) relative to preindustrial CH,. 
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Extended DataFig. 7 | Regional annual mean vertical profiles of 
instantaneous heating rate for 10xCH,sw. As in Fig. 4 from the main paper 


(which is based on the global annual mean), but broken down into the (a, b) high- 


latitude annual mean; (c, d) mid-latitude annual mean; and (e, f) low-latitude 
annual mean. Instantaneous atmospheric (a, c, e) clear-sky shortwave heating 
rate (QRS IRF,4 gray) and the corresponding clear-sky shortwave heating rate 
without water vapor shortwave absorption (QRS IRF., noH2Ov; purple) in the 
same three near-infrared bands (1.6-1.9, 2.15-2.5 and 3.1-3.85 um) that methane 


absorbs in. Instantaneous atmospheric (b, d, f) all-sky (that is, with clouds) 
shortwave heating rate (QRS IRF; gray) and the corresponding shortwave heating 
rate without water vapor shortwave absorption (QRS IRF noH2Ov; purple) in the 
samethree near-infrared bands that methane absorbs in. Also included in (a, c, e) 
isthe climatological specific humidity (SH; red) and in (b, d, f) the climatological 
cloud fraction (CLOUD; cyan). QRS IRF, QRS IRF,,, CLOUD and SH units are K day”, 
K day”, % and gkg”, respectively. 
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PDRMIP 3xCH4 Results 
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| *Error bar center is the PDRMIP model mean response 4 
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Extended Data Fig. 8 | 3xCH, PDRMIP results. Global annual mean (a) 
effective radiative forcing (ERF), (b) near-surface air temperature response, 
(c) precipitation response, and (d) apparent hydrological sensitivity to 3xCH, 
from PDRMIP experiments. The first bar in each like-colored set of three bars 
represents 4 models that include methane LW and SW radiative effects; the 
second bar represents 6 models that include the methane LW radiative effect; 
and the third bar represents the difference (that is, an indirect estimate of 
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the methane SW radiative effects). Error bars show the 1-standard deviation 
uncertainty across models (added in quadrature for the third bar). This includes 
4 models for methane LW and SW radiative effects and 6 models for methane 
LW radiative effects. Units in (a) are W m ?; units in (b) are K; units in (c) are 96; 
and units in (d) are % K”. ERF comes from PDRMIP fSST simulations; remaining 
responses are estimated from the coupled simulations. Symbols (X's) in the first 
two sets of bars for each panel represent individual PDRMIP models. 
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